Vortex lattice formation in a rotating Bose-Einstein condensate 
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We study the dynamics of vortex lattice formation of a rotating trapped Bose-Einstein condensate 
by numerically solving the two-dimensional Gross-Pitaevskii equation, and find that the condensate 
undergoes elliptic deformation, followed by unstable surface-mode excitations before forming a quan- 
tized vortex lattice. The origin of the peculiar surface-mode excitations is identified to be phase 
fluctuations at the low-density surface regime. The obtained dependence of a distortion parameter 
on time and that on the driving frequency agree with the recent experiments by Madison et al. 
[Phys. Rev. Lett. 86, 4443 (2001)]. 
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Quantized vortices have long been studied in super- 
fluid "^He as the topological defects characteristic of su- 
perfluidity 0, ||. However, the relatively high density 
and strong repulsive interaction complicate the theoreti- 
cal treatments of the Bose-Einstein condensed liquid, and 
the healing length of the atomic scale makes the exper- 
imental visualization of the quantized vortices difficult. 
The recent achievement of Bose-Einstein condensation in 
trapped alkali-metal atomic gases at ultra low tempera- 
tures has stimulated intense experimental and theoretical 
activity. The atomic Bose-Einstein condensates(BECs) 
have the weak interaction because they are dilute gases, 
thus being free of the above difficulties that superfluid 
^He is subject to. Quantized vortices in the atomic BECs 
have recently been created experimentally by Matthews 
et al. Madison et al. [|| and Abo-Shaeer et al. 

By rotating an asymmetric trapping potential, Madi- 
son et al. at ENS succeeded in forming a quantized vor- 
tex in ®^Rb BEC for a stirring frequency that exceeds a 
critical value Q . Vortex lattices were obtained for higher 
frequencies. The ENS group subsequently observed that 
vortex nucleation occurs via a dynamical instability of 
the condensate ||^. For a given modulation amplitude 
and stirring frequency, the steady state of the conden- 
sate was distorted to an elliptic cloud, stationary in the 
rotating frame, as predicted by Recati et al. |^]. An in- 
trinsic dynamical instability |p[ of the steady state trans- 
formed the elliptic state into a more axisymmetric state 
with vortices. However, the origin of that instability and 
how it leads to the formation of vortex lattices remains 
to be investigated 

The ENS group found that the minimum rotation fre- 
quency r^nuc at which one vortex appears is 0.650-"^, 
where oj^ is the transverse oscillation frequency of the 
cigar-shaped trapping potential, independent of the num- 
ber of atoms or the longitudinal frequency ojz @. There 
is a discrepancy between the observations and the the- 
oretical considerations based on the stationary solution 
of the Gross-Pitaevskii equation(GPE) 0; Vt^uc is 
significantly larger than its equilibrium estimates. 

The present paper addresses these issues by numeri- 
cally solving the GPE that governs the time evolution of 
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Here g = ATrh^a/m is the coupling constant, propor- 
tional to the ^^Rb scattering length a «5.77 nm. The 
high anisotropy of the cigar-shaped potential used in 
the ENS experiments (u;±/ujz ~ 14) may permit the 
two-dimensional analysis. We thus focus on the two- 
dimensional dynamics of Eq. (^ by assuming the trap- 
ping potential 

Vtrir) = + e,)x^ + (1 + ey)y^}, (2) 

where uj± = 2tt x 219 Hz, and the parameters ex = 0.03 
and Sy — 0.09 describe small deviations of the trap from 
the axisymmetry, corresponding to the ENS experiments 
[Q. The centrifugal term —VlLz = ifiVt{xdy — ydx) ap- 
pears in a system rotating about the z axis at a fre- 
quency f2. An important characteristic parameter of the 
two-dimensional system is C = SnNa/L, with the total 
number N of the condensate atoms and L the size of the 
system along the z axis. 

The term with 7 in Eq. introduces the dissipation. 
Although the detailed mechanism of the dissipation is yet 
to be understood, we include this term in the GPE be- 
cause of the following reasons. First, collective damped 
oscillations of the condensate have been observed by Jin 
et al. and Mewes et al. , which shows the presence of 
some dissipative mechanisms and is consistent with the 
solution of the GPE with the dissipative term ||l^ . Sec- 
ond, even if the trapping potential is rotated fast enough, 
vortex lattices will never be formed without dissipation, 
because vortex lattice correspond to local minimum of 
the total energy in the configuration space |lj] . In other 
words, the observation of the vortex lattices JM implies 
the presence of dissipation. Following Ref. [ [l3[ , we use 
7 — 0.03 throughout this work. 

In this paper, we focus on two kinds of responses of the 
condensate, namely, the response to the sudden turn on 
of the rotation of the potential, and that to the slow turn 
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FIG. 1: Time development of the condensate density after the trapping potential begins to rotate suddenly with f2 — 0.7uj±. 
The time is t =0 msec (a), 21 msec (b), 107 msec (c), 114 msec (d), 123 msec (e), and 262 msec (f). The unit for length is 
ano ~ \/h/2m(jj± = 0.512/im and the period of the trap 4.57 msec. 



on. The numerical calculations are performed with the 
implicit Crank-Nicolson method. The chemical potential 
H is adjusted at all times so as to preserve the number of 
condensate atoms. 

We first prepare an equilibrium condensate with C = 
1400 trapped in a stationary potential. Figure 1 shows 
the typical dynamics of the condensate density \tp{r, i)p 
after the potential begins to rotate suddenly with = 

0. 7uj± |l5|. The condensate is elongated along the x axis 
because of the small anisotropy of Vtr [Eq. (|^)] , and the el- 
liptic cloud oscillates. Then, the boundary surface of the 
condensate becomes unstable, exciting the surface waves, 
which propagate along the surface. The excitations are 
likely to occur on the surface whose curvature is low, i.e., 
parallel to the longer axis of the ellipse. The ripples on 
the surface develop into the vortex cores around which 
superflow circulates. Subject to the dissipative vortex 
dynamics, some vortices enter the condensate, forming a 
vortex lattice. As the vortex lattice is being formed, the 
axial symmetry of the condensate is recovered by trans- 
ferring angular momentum into quantized vortices. 

This peculiar dynamics is understood by investigating 
the phase of tp(r, t) as shown in Fig. 2. There are some 
lines where the phase changes discontinuously from black 
to white, which corresponds to the branch cuts between 
the phases and 27r. Their ends represent phase defects, 

1. e. vortices. As soon as the rotation starts, some defects 
begin to enter. When the defects are on the outskirts 
of the condensate where the amplitude |'0(r, t)\ is almost 
negligible, they neither contribute to the energy nor the 
angular momentum of the system. These defects come 
into the boundary surface of the condensate within which 
the amplitude grows up. Then the defects compete with 



each other and induce the above surface waves due to in- 
terference. There the selection of the defects starts, be- 
cause their further invasion into the condensate costs the 
energy and the angular momentum. As is well known in 
the study of rotating supcrfluid helium , the rotating 
drive pulls vortices into the rotation axis, while repulsive 
interaction tends to push them apart; this competition 
yields a vortex lattice whose vortex density depends on 
the rotation frequency. In our case, some vortices enter 
the condensate and form a lattice dependent on f2, while 
excessive vortices are repelled and escape to the outside. 
Remarkably, the phase profile of Fig. 2(b) reveals that 
the repelled vortices also form a lattice on the outskirts of 
the condensate. Since they cannot be seen in the corre- 
sponding density profile of Fig. 1(f), they may be called 
"ghost" vortices. 



(a) (b) 
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FIG. 2: Phase profile of tp: (a) and (b) corresponds to Fig. 
1 (c) and 1 (f). The value of the phase varies continuously 
from (black) to 2tt (white). The unit for length is the same 
as that of Fig. 1 
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The distortion of the condensate to an eUiptic cloud 
was theoretically studied by Recati et al. 0. Assuming 
the quadrupolar velocity field v(r) = aSI{xy), they ob- 
tained the distortion parameter a = il{R^ — Ry)/{Rl + 
Ry) in the steady states as a function of fi, where R^ 
and Ry are the sizes of the condensate along each direc- 
tion. Madison et al. observed that, after the rotation of 
n = 0.7wj_ starts suddenly, a oscillates during a few hun- 
dred milliseconds and then falls abruptly to a value below 
0.1 when vortices enter the condensate from its boundary 
surface Figure 3 shows the oscillation of a and the 
increase of the angular momentum £z/h per atom in our 
dynamics of Fig. 1. This figure closely resembles Fig. 3 
of Rcf. and our scenario is thus consistent with the 
experimental results. 



The dispersion relation of the surface waves |16 of a 
rotating condensate is obtained as follows. The substitu- 
tion of the Madelung transformation ip — j-^le** into Eq. 
(§ yields 



dt 



+ V- [|7/'|2(v-f] X r)] ^0, (3a 



dv 
~dt 



m— = -V{6fi+^{yr-nxrf-'^n\^}, (3b) 



where v = {h/m)V0 and = Vtr+g\ip\^- ^T^V^jV-l - 
fi. By linearizing Eq. (3), we obtain the dispersion re- 



i?([j| -n^)k with the 



lation of the surface waves cu'^ = i?(ij-'x 
condensate size R. The dependence of on i? shows 
that the surface with lower curvature is excited more eas- 
ily, which is clearly seen in Fig. 1. The group velocity 

agrees with the propaga- 
tion velocity of the surface wave in our simulation. As 
discussed by Al Khawaja et al. the surface waves 

are connected with the low energy excitations studied 
by Stringari and Dalfovo et al. |]l7| , and Isoshima and 
Machida 

We have examined the critical frequency at which vor- 
tices can enter with C — 420 corresponding to the ex- 
perimental condition and found that only when il 
is larger than Qci — 0.57a;j^ can vortices enter the con- 
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FIG. 3: Time evolution of the distortion parameter a (solid 
line) and the angular momentum £z/h per atom (dashed line) 
corresponding to the dynamics of Fig. 1. 



densate and form a lattice. This critical frequency is 
closer to the observed value ilnnc — 0.65w_l, than the 
values obtained in previous literature. The critical fre- 
quencies have been studied from thermodynamic or sta- 
bility arguments |l^ , but they are generally much 
smaller than the observed value. It should be noted that 
these critical frequencies only give the necessary condi- 
tion which enables a vortex to exist stably at the center 
of the trap. Actually, vortices should be nucleated at the 
boundaries and come into the condensate; the condition 
of ri, which realizes such nonlinear dynamics, may be 
generally different from that obtained from the stability 
arguments. Isoshima and Machida examined the local in- 
stability of nonvortex states toward vortex states within 
the Bogoliubov theory jl^. They note that this local 
instability could give higher critical frequencies than the 
global instability, which comes from comparing their to- 
tal energy, for example, about 0.6a;^ for the ENS exper- 
imental condition |l9|]. The stability analysis cannot an- 
swer what happens when the system goes over the linear 
region. Our paper reveals the importance of the nonlin- 
ear dynamics beyond the stability analysis. The detailed 
studies of the dependence of the number of vortices on 
that of r^ci on the number of condensate atoms pC| ] 
will be reported elsewhere. 

Finally, we study the evolution of the condensate in a 
potential rotating with a time-dependent frequency Q{t). 
The experiments were made by Madison et al. H. The 
condensate at every moment was in a stationary state 
with an elliptic cloud obtained by Recati et al. under 
the Thomas- Fermi limit Q, since the time dependence 
was slow enough. One of the important observations 
was that, for an ascending ramp, the condensate follows 
branch I of the nonvortex state until Qc — 0.75a;_L- Be- 
yond this critical frequency, the condensate changes from 
branch I to the more axisymmetric state, nucleating vor- 
tices. This critical frequency agrees with the frequency at 
which branch I becomes dynamically unstable how- 
ever, the detailed nature of the bifurcation has remained 
to be clarified. We study the response of the conden- 
sate of C = 420 for the ascending frequency ^l{t) = ct 
with c = 2tt X 400 Hz/s; they are about the same values 
as those used in the ENS experiment. The obtained re- 
sults are consistent with the experimental ones, as shown 
in Fig. 4. As increases, the condensate becomes 
gradually elliptic following branch I, when the nucleated 
"ghost" vortices come to the surface of the condensate. 
When is raised to the critical value ilc2 — 0.75a;_L, 
the vortices jostle with each other and excite the surface 
waves, reflecting the dynamical instability. Then, sim- 
ilar to the first set of calculations, some vortices enter 
the condensate, form a lattice and increase the angular 
momentum, while the condensate recovers the axisym- 
metry apart from branch I. The small discrepancy of the 
trace of our calculation from branch I in Fig. 4 may be 
attributed to a deviation from the Thomas- Fermi limit; 
the similar discrepancy is reported by the Madison et al. 
too 0. 
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FIG. 4: Development of a for the ascending frequency Q,{t). 
The sohd curve describes our numerical result, while the 
dashed curve represents the stationary nonvortex state ob- 
tained by Recati et al, See the text for more details. 

It is necessary to compare two critical frequencies 
ilci — 0.57w^ and 17^2 — Q-lhoj^. What causes the differ- 
ence? In the first set of calculations where the potential 
begins to rotate abruptly with a frequency fJ, the conden- 
sate first deviates much from the stationary state a{^) 
with the quadrupolar velocity field. As shown in Fig. 3, 
it seemingly relaxes to that state through the damped 
oscillation of the large amplitude, but it changes on the 
way to another state with a vortex lattice which has lower 
energy than the nonvortex state a{Vl). The critical fre- 
quency f2ci means that the state with vortices can be 
realized dynamically from the nonvortex state only when 

> r2ci • In the second calculation, as f2 increases slowly, 
the condensate follows quasistatically the nonvortex state 
a{Vl) of branch I, even after it becomes metastable for 

> r^ci- Only when exceeds f2c2, the condensate can 
change from this state to the vortex state, reflecting the 
dynamical instability 

Recently Garcia-RipoU et al. and Feder et al. study 
this system by three-dimensional analysis of the GPE 



JlSf , taking account of the bending of vortex lines in 
the cigar-shaped potential. A bent vortex moves with 
the self-induced velocity proportional to its curvature 
0, which is lacking in the two-dimensional dynamics. 
However, due to the high anisotropy of the cigar-shaped 
potential, vortices are aligned almost straight along the 
symmetry axis. It is therefore unlikely that the bend- 
ing effect alters the results of our work qualitatively and 
significantly. 

In conclusion, we have studied the dynamics of vortex 
lattice formation of a rotating trapped BEG by numeri- 
cally solving the two-dimensional GPE, and obtained the 
following physical picture. When the trapping potential 
begins to rotate at a sufficiently fast frequency, the con- 
densate is distorted to an elliptic shape and oscillates. 
Then, its boundary surface becomes unstable, exciting 
surface waves. The origin of these ripples is identified 
to be the interference of violent phase fluctuations that 
occur on the outskirts of the condensate; and some of the 
surface ripples develop into the vortex cores which then 
form a vortex lattice. The critical frequency for this pro- 
cess is found to be fici — 0.57w^ for the condition at the 
ENS group. On the other hand, if the rotation frequency 
is raised slowly from zero, the condensate becomes el- 
liptic and follows the nonvortex elliptic state with the 
quadrupolar flow 0. When exceeds ilc2 — 0.75a;^, 
the condensate deviates from the state, nucleating vor- 
tices via the dynamical instability Q . The whole dynam- 
ics is concerned with the "ghost" vortices, which are the 
phase defects outside the condensate. The obtained re- 
sults are consistent with the experimental results of the 
ENS group |,|. 
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